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ABSTRACT 

We report observations of the hot post-asymptotic giant branch star ZNG 1 in the globular cluster 
M5 (NGC 5904) with the Far Ultraviolet Spectroscopic Explorer (FUSE). From the resulting spectrum, 
we derive an effective temperature Toff = 44300 ± 300 K, a surface gravity log g = 4.3 ±0.1, a rotational 
velocity w sin i = 170±20 km s"^, and a luminosity log L/L© = 3.52±0.04. The atmosphere is helium-rich 
(Y = 0.93), with enhanced carbon (2.6% by mass), nitrogen (0.51%) and oxygen (0.37%) abundances. 
The spectrum shows evidence for a wind with terminal velocity near 1000 km and an expanding shell 
of carbon- and nitrogen-rich material around the star. The abundance pattern of ZNG 1 is suggestive 
of the "born-again" scenario, whereby a star on the white-dwarf cooling curve undergoes a very late 
shell flash and returns to the AGB, but the star's rapid rotation is more easily explained by a previous 
interaction with a binary companion. 

Subject headings: globular clusters: individual (NGC 5904) — stars: AGB and post-AGB — stars: 
individual (CI* NGC 5904 ZNG 1) — stars: mass-loss — ultraviolet: stars 



1. INTRODUCTION 

Post-asymptotic giant branch (PAGB) stars represent a 
brief (lO'^-lO^ yr) phase of stellar evolution during which 
stars move rapidly across the HR diagram at constant lu- 
minosity from the cool tip of the AGB to the beginning of 
the hot white-dwarf cooling curve. As such, they are keys 
to both the past and the future: their atmospheric abun- 
dances provide important information about the mixing 
and mass-loss processes at work in AGB stars, while their 
masses yield constraints on the masses of white dwarfs cur- 
rently forming in globular clusters. To investigate these 
issues, we observed the PAGB star ZNG 1 (Table 1) in the 
globular cluster M5 (NGC 5904) with the Far Ultraviolet 
Spectroscopic Explorer (FUSE). 

First identified as a UV-bright star by Zinn, Newell, & 
Gibson (1972) on the basis of its extreme U — V color, 
ZNG 1 was observed with lUE by Bohlin et al. (1983), 
who suggested that an apparent N IV A1487 emission fea- 
ture is the signature of a planetary nebula (PN) . Combin- 
ing this spectrum with other archival lUE data, de Boer 
(1985) argued that the apparent N IV feature is spuri- 
ous and pointed out the P-Cygni shape of the strong N V 
AA1238, 1242 doublet. To test the PN hypothesis, Napi- 
wotzki & Heber (1997) used the Hubble Space Telescope to 
obtain a WFPC2 Hq image and GHRS ultraviolet spectra 
of the star. The Ha image shows a bright G-type star only 

1 Based on observations made with the NASA-CNES-CSA Far Ultr 
Hopkins University under NASA contract NAS5-32985. 



0'.'5 from ZNG 1, but no evidence for extended emission 
around the star. The GHRS spectra clearly show P-Cygni 
profiles in the N V doublet and discrete absorption compo- 
nents (DACs), blueshifted by - 900 km s'^, in both N V 
and C IV AA1548, 1550. No nebular emission is present. 

2. OBSERVATIONS AND DATA REDUCTION 

FUSE consists of four separate optical systems. Two 
employ LiF optical coatings and are sensitive to wave- 
lengths from 990 to 1187 A, while the other two use SiC 
coatings, which provide reflectivity to wavelengths as short 
as 905 A. The four channels overlap between 990 and 1070 
A. The FUSE flux calibration, based on theoretical mod- 
els of white-dwarf stellar atmospheres, is believed accurate 
to about 10%. For a complete description of FUSE, see 
Moos et al. (2000) and Sahnow et al. (2000). 

The FUSE spectrum of ZNG 1 in M5 (data set 
A1080303) was obtained in 4 separate exposures on 2000 
July 15. The total integration time was 3681 s, all of it dur- 
ing orbital night. All observations were made through the 
30" X 30" (LWRS) aperture. Archival HST/WFPC2 FUV 
images of the cluster core confirm that no other UV-bright 
stars fell within the spectrograph aperture. The data 
were reduced using version 2.3 of the CalFUSE calibration 
software pipeline. For each FUSE channel, the spectra 
from individual exposures were cross-correlated, shifted to 

t Spectroscopic Explorer. FUSE is operated for NASA by the Johns 
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a common wavelength scale, weighted by their exposure 
time, and averaged. The spectrum from each channel was 
then cross-correlated with a synthetic molecular-hydrogen 
spectrum and shifted to match. The wavelength scale de- 
fined in this way is at rest with respect to the interstellar 
medium (ISM), which we take as a proxy for the local 
standard of rest. Finally, the spectra were binned by four 
detector pixels, or about 0.025 A, approximately half of 
the instrument resolution. 

3. SPECTRAL ANALYSIS 

3.1. Synthetic Spectra and Model Fitting 

Model atmospheres are computed using the program 
TLUSTY (Hiibeny & Lanz 1995), which employs a sophis- 
ticated treatment of non-LTE (local thermodynamic equi- 
librium) line blanketing to account for the opacity sources 
important in stellar atmospheres. In our models, the line 
and continuum transitions of H, He, C, and N (a total of 
10 ions and 220 energy levels) arc treated in non-LTE with 
detailed photoionization cross sections that include reso- 
nances; other species, assumed to liavc scaled-solar abun- 
dances with [Fe/H] equal to the cluster mean (Table 1), 
are treated in LTE. All models assume a microturbulent 
velocity of 2 km s"'^. Synthetic spectra are generated us- 
ing the program SYNSPEC and rotationally broadened 
using the program R0TIN3, which assumes constant limb 
darkening. The application of these codes to another sdO 
star, BD -|-75°325, is described in Lanz, Hubeny, & Heap 
(1997). 

Model spectra are fit to the data using the nonlinear 
curve-fitting program SPECFIT (Kriss 1994), which per- 
forms a minimization of the model parameters. Red- 
dening is modeled with a Cardelli, Clayton, & Mathis 
(1989) extinction curve, extrapolated to the Lyman limit, 
assuming E{B-V) = 0.03 (Table 1) and Ry = 3.1. Er- 
ror bars for a particular parameter are derived by fixing 
that parameter at the best-fit value, then raising it. while 
allowing the other model parameters to vary freely, until 
rises by 1, which corresponds to a 1 cr deviation for a 
single interesting parameter (Avni 1976). The error bars 
quoted in Table 2 are purely statistical and do not include 
systematic errors. 

3.2. Atmospheric Parameters 

The most prominant feature in the FUSE spectrum of 
ZNG 1 is the C III* multiplet at 1175 A. The half dozen 
individual lines of this feature are generally well resolved 
in FUSE spectra, but in ZNG 1 arc blended into a broad 
absorption trough. As shown in Fig. 1, this feature is well 
reproduced by a rotational velocity vsini = 170 km s"-'^, a 
carbon mass fraction of 3%, an effective temperature Teg 
= 45 kK, and a radial velocity Vlsr = +60 km s~^, a value 
consistent with the cluster velocity (Table 1). 

We determine the star's surface gravity and helium 
abundance from simultaneous fits to its helium and hydro- 
gen lines. As shown in Fig. 2, every second He II Balmer 
line is paired with an H I Lyman line, but lies ~ 0.5 A 
to the blue. The lines are unresolved because of the star's 
high rotational velocity, but our models show that these H- 
He blends shift to the blue as the helium abundance rises 
and the hydrogen abundance falls. (Note that the stel- 
lar features are redshifted relative to the ISM.) Though 



the cores of all three features are contaminated — the H- 
He blends by interstellar H I and the He II A933 line by 
S VI from the stellar wind — their wings can constrain both 
the stellar surface gravity and the helium abundance. 

We model the region between 930 and 940 A, exclud- 
ing the cores of all three stellar features (Fig. 2). (Model 
fits to the ISM lines confirm that all regions of significant 
foreground H I absorption are excluded.) The remaining 
interstellar lines are fit with simple Gaussians. We gen- 
erate a grid of models with surface gravities log g = 4.0, 
4.5, and 5.0 and helium fractions from 10 to 90% (in steps 
of 10%) by number. Within SPECFIT, we linearly inter- 
polate between models to produce a fine grid with steps 
of 0.1 dex in logg. Assuming Tofr = 45 kK and Vlsr = 
+60 km s^^, wo find a surface gravity of log g = 4.3 ±0.1 
and a helium fraction of 85 ± 5% by number. Repeating 
this test with T^g = 46 kK models yields best-fit values 
of log g = 4.4 ±0.1 and a helium fraction of 85 ± 5%, so 
wc need not increase our error bars to account for small 
uncertainties in the effective temperature. Given the low 
hydrogen abundance of this star, we quote all abundances 
as mass fractions; that of helium is 93 ± 2% (Table 2). 

Adopting these stellar parameters, we fit model spectra 
to four photospheric nitrogen features, N IV A955, N III 
A980, N III A990 (red wing), and N III A992, and derive 
an effective temperature of Tes = 44300 ± 300 K and a 
nitrogen mass fraction of 0.51 ± 0.05%. Fits to three of 
these lines arc presented in Fig. 1. The giants in M5 show 
a wide range of nitrogen abundances, the highest being 
[N/Fe] = +1.38 (Cohen, Briley, & Stetson 2002), but the 
nitrogen abundance of ZNG 1 is 4 times greater. Return- 
ing to the C III* multiplet, we obtain a rotational velocity 
vsini = 170 ± 20 km s^^ and a carbon mass fraction of 
2.58 ± 0.20%, more than 160 times greater than the most 
carbon-rich cluster giant studied by Cohen et al. (2002). 

The star's FUSE spectrum shows no strong oxygen fea- 
tures, so we fit the wavelength region between 1090 and 
1104 A, where the continuum is fairly flat and where nu- 
merous weak O III and O IV lines combine to form two 
broad absorption troughs centered at 1093 and 1099 A. 
We derive an oxygen mass fraction of 0.37 ±0.32%. While 
the oxygen abundance is the least well-determined of the 
CNO abundances in ZNG 1, it is certainly larger than 
the maximum value observed in M5 giants, where [0/Fe] 
ranges from -0.55 to +0.47 (Ivans et al. 2001). The star 
also appears to be enhanced in S, Si, and P; we will deter- 
mine their abundances in a future paper. 

3.3. Wind and Shell 

The FUSE spectrum of ZNG 1 shows broad P-Cygni 
profiles and blueshifted discrete absorption components 
(DACs) in both components of the O VI doublet. The 
O VI profiles are well fit by models with a terminal ve- 
locity of Voo = 1000 km s~^. The DACs, which are gen- 
erally attributed to density enhancements in the stellar 
wind, are blueshifted by 960 km s~^ relative to the star, 
consistent with the DACs seen by GHRS (Napiwotzki & 
Heber 1997). The spectrum also shows narrow absorption 
features of H I, C III A977.02, C II A1036.33, and N II 
A1083.99 blueshifted by ~ 120 km s'^ relative to the ISM 
and ~ 180 km s~^ relative to the star (two blue-shifted 
H I features are marked in Fig. 2; see also Fig. 1 of Dixon, 
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Brown, & Landsman 2003). These features may represent 
a shell of material expanding about the star. A complete 
analysis of the star's wind and shell features is forthcom- 
ing. 

3.4. Stellar Parameters 

Using the cluster distance and the ratio of observed to 
synthetic spectra (Dixon, Davidsen, & Ferguson 1994), 
we derive a stellar luminosity logL/LQ = 3.52 ± 0.04. 
The star lies on the Schonberner (1983) PAGE evolution- 
ary tracks and should have a mass M^, ^ 0.565. The 
mass calculated from the parameters of our best-fit model, 
M* = 0.69 ± 0.17 Mq, is consistent with the Schonberner 
prediction. 

4. DISCUSSION 

Hot PAGE stars in globular clusters differ from their 
field counterparts in that they rarely show evidence for 
third dredge-up (by which the products of helium burning 
are brought to the surface) and are seldom accompanied by 
PNe (Jacoby et al. 1997). These characteristics can be un- 
derstood qualitatively as a result of the expected low core 
mass (< O.55M0) of globular cluster AGB stars. While the 
minimum core mass for third dredge-up is still a subject of 
active research (Karakas, Lattanzio, & Pols 2002), nearly 
all current AGE models agree that third dredge-up should 
not occur in stars with initial mass ^0.8 Mq. The rar- 
ity of PNe can be explained by the lower luminosities and 
longer crossing times of PAGE stars with low core masses; 
the nebula disperses before the central star becomes hot 
enough to ionize it (Jacoby et al. 1997). Alves, Eond, 
& Livio (2000) therefore suggest that binary evolution is 
necessary for the formation of PNe in globular clusters. 
Interestingly, the best-studied PN in a globular cluster, 
K648 in M15, shows clear evidence for third dredge- up 
in its central star (Rauch, Heber, & Werner 2002), but 
not the strong helium and nitrogen enhancements seen in 
ZNG 1. 

One way to produce a hydrogen-deficient star with en- 
hanced CNO abundances is for the star to experience a 
very late shell flash while descending the white dwarf cool- 
ing curve and subsequently return to the AGB (Herwig 
et al. 1999). This "born-again" scenario might explain the 
absence of a PN surrounding ZNG 1, because the star is so 
old (since it is on a second crossing), that any nebula has 
long since dispersed. For the born-again scenario to work 
for ZNG 1, the hydrogen must not be completely burned, 
and the carbon enhancement must not be too high. Lawlor 
& MacDonald (2003) provide a born-again model of an 
initial 1 Mq star with Z = 0.001 that predicts stellar pa- 
rameters similar to those of ZNG 1 (see their Table 3), 
although the CNO abundances are about a factor of three 
larger. The difiiculty with the born-again scenario (or any 
single-star evolution scenario) is providing an explanation 
for the observed rapid rotation. One possibility: to ex- 
plain the observed rotational velocity of up to 40 km s~^ 
in some globular horizontal-branch stars (Eehr 2003), Sills 
& Pinsonneault (2000) postulate that red giants may have 
a rapidly rotating core. If so, then perhaps its angular 
momentum can be transferred to the surface during the 
convective mixing in the final shell fiash. 

Might the properties of ZNG 1 be explained through 



interaction with another star? First, we note that its 
very blue B — V color indicates that any current compan- 
ion must be quite faint. Given the proximity of the star 
to the projected cluster center (4'.'2) and the moderately 
high stellar interaction rate in M5 (Pooley et al. 2003), 
one possiblity is a stellar collision, which is thought to be 
the primary production mechanism for the blue stragglers 
found at the centers of globular clusters. CoUisionally- 
created blue stragglers are expected to be born with a 
high angular momentum. Shara, Saffer, & Livio (1997) 
measure a rotation rate of 155 ± 55 km s~^ and estimate 
a mass of 1.7 ± 0.4 Mq for a blue straggler in 47 Tuc. A 
star of this mass would likely experience third dredge-up 
on the AGB, but would also be expected to eject a PN, 
whic;h is not secin for ZNG 1. Furthermore, it is unclear 
how much of a blue straggler's angular momentum would 
survive its post-main sequence evolution. 

A second scenario, in which a low-mass companion spi- 
rals into the AGB progenitor, has been proposed to explain 
K648 in M15 by Alves et al. (2000) and discussed in gen- 
eral by De Marco & Soker (2002). The increased AGE 
mass due to the ingestion of the companion might be suf- 
ficient to allow third dredge-up to occur, but additional 
mechanisms such as shear or rotationally-induced mixing 
are required to create a hydrogen-deficient star (De Marco 
& Soker 2002). A merger process is thought to be the 
origin of the high rotational velocity of FK Comae stars 
(usini = 162.5 ± 3.5 km for the prototype FK Co- 
mae; Huenemoerder et al. 1993). More detailed analysis 
is required to assess the viability of a merger scenario for 
ZNG 1; whether the merger remnant should eject a PN is 
of particular interest. 

Pauldrach et al. (1988) compute terminal velocities and 
mass-loss rates for the winds of (Pop. I) central stars of 
PNe. Assuming an evolutionary mass of 0.565 Mq (§3.4) 
and temperature of Toff = 45 K, their models predict mass- 
loss rates of a few times 10~^ Mq yr~^. Given this modest 
wind, it appears unlikely that ZNG 1 will dissipate sig- 
nificant angular momentum in the brief time (~ 15,000 
years, according to the 0.565 Mq track of Schonberner 
1983) prior to its arrival on the white-dwarf cooling curve. 
Previous studies have yet to find white dwarfs with such 
rapid rotation (see review by Kawaler 2003) but the ex- 
ample of ZNG 1 suggests that further searches may be 
fruitful. 

This research has made use of NASA's Astrophysics 
Data System Bibliographic Services and the Catalogue 
Service of the CDS, Strasbourg, France. We thank \. 
Hubeny for supplying his spectral synthesis codes and A. 
Sweigart for illuminating discussions. We acknowledge the 
outstanding efforts of the FUSE P.L team to make this 
mission successful. This work is supported by NASA grant 
NAG 5-10405. 
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Fig. 1. — Nitrogen and carbon features in the FUSE spectrum of ZNG 1 in M5. Data are binned to 0.025 A and plotted as histograms. 
Error bars and synthetic spectra are overplotted. Top panel: The N IV A955, N III A990, and N III A992 lines overplotted by models with 
Teff = 40 (dotted line), 45 (solid line), and 50 kK (dashed line). The narrow absorption features are due to interstellar H2 and O I. Bottom 
panel: The C III* A1175 multiplet overplotted by models with rotational velocity vsini = 10 (gray line), 100 (dotted line), 200 (solid line), 
and 300 (dashed line) km s~^. 
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Fig. 2. — Hydrogen and helium lines in the FUSE spectrum of ZNG 1 in M5. Data are taken from the two SiC detector segments, binned 
to 0.025 A, and plotted as a histogram. Error bars and synthetic spectra are overplotted. The thick, solid Uncs below the spectrum mark 
the spectral regions included in the model fits. The models have T^g = 45 kK, logg = 4.5, and helium fractions of 10% (dotted line), 40% 
(dashed line), and 90% (solid line) by number. As the helium fraction increases, the helium- hydrogen blends at 931 and 938 A shift to the 
blue and the He II A934 line becomes stronger. Above the spectrum are marked two pairs of foreground H I absorption lines; the bluer pair 
may represent a shell of material expanding about the star. 
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Table 1 
Observational Parameters 



Parameter 


Value 


Reference 


Spectral type 


sdO 


1 


V 


14.54 


2 


B-V 


-0.32 


2 


E{B-V) 


0.03 


3 


Distance (kpc) 


7.5 


3 


[Fc/H] cluster 


-1.27 


3 


VLSli.cliL.Mcr 


()."). 7 


;] 



References. — (1) Napiwotzki & 
Hebcr 1997, (2) Piotto et al. 2002, 
(3) Harris 1996 (on-line version dated 
February 2003). 
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Table 2 



Derived Stellar Parameters 



Parameter 


Stellar Value 




44300 ± 300 K 


logs' 


4.3 ±0.1 


vsini 


170 ± 20 km s-i 


^LSR 


+60 ±4 km s-i 


logL/L© 


3.52 ± 0.04 


M/Mq 


0.69 ±0.17 


Y 


0.93 ±0.02 


C abundance 


2.58 ±0.20% 


N abundance 


0.51 ±0.05% 


abundance 


0.37 ±0.32% 



Note. Abundances are quoted 
as mass fractions. 



